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Abstract. In situ measurements of hydrogen, nitrogen, and 
chlorine radicals obtained in the lower stratosphere during SPADE 
are compared to results from a photochemical model that 
assimilates measurements of radical precursors and environmental 
conditions. Models allowing for heterogeneous hydrolysis of 
N 2 0 5 agree well with measured concentrations of NO and CIO, 
but concentrations of H0 2 and OH are underestimated by 10 to 
25%, concentrations of N0 2 are overestimated by 10 to 30%, and 
concentrations of HC1 are overestimated by a factor of 2. 
Discrepancies for [OH] and [H0 2 ] are reduced if we allow for 
higher yields of 0( l D) from 0 3 photolysis and for heterogeneous 
production of HN0 2 . The data suggest more efficient catalytic 
removal of O3 by hydrogen and halogen radicals relative to 
nitrogen oxide radicals than predicted by models using 
recommended rates and cross sections. Increases in [0 3 ] in the 
lower stratosphere may be larger in response to inputs of NO, 
from supersonic aircraft than estimated by current assessment 
models. 

Introduction 

Data collected during the Stratospheric Photochemistry, Aero- 
sols, and Dynamics Expedition (SPADE) in May 1993 offer a 
unique opportunity to investigate the chemistry of HO„ NO x , and 
CIO, radicals in the lower stratosphere. We assimilate into a pho- 
tochemical model data for pressure, temperature, aerosol surface 
area, long-lived precursors of free radicals (0 3 , H 2 0, CH 4 , NO,, 
and Cl, inferred from CC1 3 F and C0 2 ), total column ozone, and 
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planetary albedo obtained by the ER-2 aircraft on latitude tran- 
sects from 15 to 58°N. Concentrations of radicals (OH, H0 2 , NO, 
N0 2 , CIO) and HC1 are computed and compared to observations 
to test our understanding of key photochemical processes. 

Theory and Observation 

We focus on segments of two flights, to Moffett Field, CA 
(37.4°N) from the south on 14 May 1993 (930514) and from the 
north on 18 May 1993 (930518), that provided the most complete 
set of measurements during SPADE. Each segment began with a 
descent from -30 mb to -120 mb near noon, followed by a climb 
and cruise at maximum altitude (-50 mb), and a second descent 
and ascent near the midpoint. Fig. 1 shows measurements (used 
as inputs to the photochemical model) of pressure, temperature, 
aerosol surface area, and long-lived radical precursors. 

The aerosol surface area was measured by the Focused Cavity 
Aerosol Spectrometer (FCAS) and the Forward Scattering Spec- 
trometer Probe (FSSP), optical particle counters that determine the 
size of ambient particles with diameters between 0.1 and 0.9 pm, 
and diameters between 0.4 and 20 pm, respectively. The surface 
area estimated by the FSSP for the overlap region between the two 
instruments was approximately 50% largeT than the area measured 
by the FCAS. We have adopted an estimate of aerosol surface 
area for the size range 0.1 to 20 pm using data from the FCAS for 
the overlap region. Model results are relatively insensitive to this 
assumption, since total surface area changes by -20% depending 
on which measurement is used for the overlap region. 

Measurements of [NO,] were not available for the latter por- 
tion of 930514 and were estimated from data for [N 2 0] using the 
compact correlation observed between these species. Concentra- 
tions of inorganic chlorine (Cl,) were computed from in situ meas- 
urements of halogenated source gases and [C0 2 ] [Woodbridge et 
al., 1994]. We adopted the mean of H 2 0 concentrations reported 
by two instruments since the -15% systematic difference between 
the measurements is not understood. Other inputs to the model 
include measured [CH*], [CO] estimated from its correlation with 
[0 3 ] [Murphy et al., 1993], and [C 2 H*] and concentrations of inor- 
ganic Br estimated from their correlations with [N z O] [D. Blake, 
private communication, 1994; Schauffler et al., 1993]. 

Concentrations of HO, ([OH]+[H0 2 ]) are sensitive to the 
column abundance of Oj due to its role in regulating production of 
0( J D); partitioning of NO, ([NO]+[N0 2 ]) is sensitive to planetary 
albedo due to photolysis of N0 2 . Photolysis rates were calculated 
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LATITUDE (N) 

Figure I. Values of pressure, temperature, aerosol surface area, 
and mixing ratios of 03 . NO,. Cl,, and H 2 0 obtained during the 
return segment of ER>2 flights 9305 14 and 930518. Dashed lines 
indicate values inferred from other tracers (see text). 


using a radiative transfer model that accounts for Rayleigh and 
aerosol scattering [Prather, 1981] using ultraviolet albedo 
specified from observations by the Total Ozone Mapping Spec- 
trometer (TOMS). Profiles for O 3 were scaled to the column 
abundances measured by TOMS assuming changes occur only for 
altitudes < 24 km (26 km for latitude > 50°N). Profiles for tem- 
perature and aerosol extinction for May 1993 were obtained from 
the National Meteorological Center and from observations by the 
Stratospheric Aerosol and Gas Experiment n, respectively. 

The photochemical model includes 35 reactive species and 
approximately 220 chemical reactions. Reaction rates and absorp- 
tion cross sections were adopted from DeMore et al. [1992], 
unless noted otherwise, except for extrapolations to longer 
wavelengths of cross sections for H 2 0 2 , HN0 2 , and HN0 4 . A 
reaction probability of 0.1 was used for the heterogeneous hydro- 
lysis of N 2 0 5 [DeMore et al., 1992]. The formulation of Hanson 
et al. [1994] was adopted for the rate of the heterogeneous hydro- 
lysis of CINO 3 , which has a negligible effect on model results for 
temperatures encountered during these flights. 

A Newton-Raphson iterative technique was used to find the set 
of diumally varying species concentrations that repeats every 24 


hours [Logan et al., 1978]; i.e., we assumed that each species 
reached a balance between production and loss, over 24 hours, for 
the temperature, pressure, and latitude at which an air parcel was 
sampled. Time constants for equilibration are «1 day for HO x , 
~1 week for NO,, and -1 month for HC1. Winds were weak dur- 
ing May, the flow was primarily zonal, and the atmosphere does 
not have strong latitudinal gradients in the [CINO 3 ] 10 [HC1] ratio; 
therefore, we expect on average relatively small deviations from 
steady state for [NO x ] and [ClNOj], 

Model results appropriate for the solar zenith angle at the time 
of measurement are compared to observations in Fig. 2. Calcula- 
tions neglecting heterogeneous chemistry underestimate concen- 
trations of OH, H0 2 , and CIO and overestimate concentrations of 
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Figure 2. Measurements (points) of [OH], [HO^, [NO], [N0 2 ], 
[CIO], and [HQ] for 930514 and 930518 compared to model 
simulations using rates and cross sections from DeMore et al. 
[1992]: neglecting all heterogeneous processes (. JPL-Gas , blue 
dotted line); including heterogeneous hydrolysis of N 2 O s and 
CINO 3 ( JPL-Het , blue solid line); and same as JPL-Het except 
allowing for updated absorption cross sections for HN0 3 and 
quantum yields of 0( 1 D), and the heterogeneous production of 
HN0 2 (Model B; red line). Measurements of [CIO] and [NO] 
(obtained every 16 and 1 seconds, respectively) have been 
smoothed using a 1 minute median filter. 
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NO and N0 2 by nearly a factor of 2, significantly larger than the 
uncertainty of measurement for each gas (2c measurement uncer- 
tainties, including random and systematic effects, are estimated to 
be 30% for [OH], 40% for [H0 2 ], 20% for [NO], 25% for [N0 2 ], 
25% for [CIO], and 20% for [HC1]). 

Calculated concentrations of radicals are within the uncer- 
tainty of measurement for OH, H0 2 , NO, N0 2 , and CIO if we 
allow for heterogeneous hydrolysis of N 2 0 5 and C1N0 3 ( JPL-Het 
in Fig. 2). Hence simultaneous observations of HO,, NO,, and 
CIO, radicals demonstrate the importance of heterogeneous reac- 
tions on sulfate aerosols (primarily hydrolysis of N 2 O s ) in regulat- 
ing the abundance of radicals in the lower stratosphere, consistent 
with previous studies using in situ measurements of [CIO] [King 
et al., 1991] and [NO] [Fahey et al., 1993]. 

Values of [OH] and [H0 2 ] observed during early morning and 
late afternoon at 37°N are represented poorly by the JPL-Het 
model (discrepancies exceed a factor of 3 near sunrise) [Wennberg 
et al., 1994; Salawitch et al., 1994]. Concentrations of HO, during 
mid-day are underestimated by -20%. Concentrations of N0 2 are 
overestimated by -25% and [NO] is slightly underestimated, sug- 
gesting possible errors in the partitioning of NO, [Jaegl6 et al., 
1994], Differences between observed and calculated concentra- 
tions of HO2 and NO 2 during mid-day suggest contributions from 
HO, and halogen cycles to the removal rate of 0 3 in the lower 
stratosphere are larger than predicted by models that use recom- 
mended rates and cross sections, although these differences lie 
within the 2o uncertainty of the measurements. 

Wennberg et al. [1994] proposed heterogeneous decomposi- 
tion of HNO4 as a source of HN0 2 to account for the rapid rise in 
[OH] and [H0 2 ] at sunrise. This reaction would increase [HO,] at 
mid-day due to reduction of the OH+HNO4 sink. Michelsen et al. 
[1994] suggested that photolysis of vibrationally and rotationally 
excited states of 0 3 result in quantum yields for 0(*D) approxi- 
mately 25% higher than values recommended by DeMore et al. 
[1992], leading to enhanced production of HO, from the reaction 
of O^DFtT^O. Recent measurements of the temperature- 
dependent absorption cross section for HN0 3 reported by Burk- 
holder et al. [1993] result in a 20% reduction in the photolysis rate 
of HNO3 for a temperature of 220 K compared to models using 
cross sections from DeMore et al. [1992], lowering calculated 
[NO x ] and [HO x ], and raising [CIO]. 

Agreement between theory and observation for [OH] and 
[H0 2 ] is improved if we adopt absorption cross sections for HN0 3 
from Burkholder et al. [1993] and quantum yields for 0( 1 D) from 
Michelsen et al. [1994], and allow for the heterogeneous decom- 
position of HN0 4 to form HN0 2 with a reaction probability of 0.2 
(A iode l B % Fig. 2). This model results in an excellent simulation of 
the diurnal variation of [OH] and [H0 2 ] observed at 37°N 
[Salawitch et al., 1994]. 

Comparisons between observed and calculated [NO,] and 
[CIO] are unable to distinguish between various heterogeneous 
models. Calculated values of [NO] and [N0 2 ] are sensitive to 
uncertainties in the rates of processes that regulate the ratio of NO 
to NO2 (including the photolysis rate of N0 2 ) and the ratio of NO, 
to NO, (including aerosol surface area, the reaction probability of 
hydrolysis of N 2 0 3 , and the influence of variations in latitude of 
air parcels during the 2 weeks prior to measurement) [Kawa et al., 
1993; R. Cohen et al., in preparation, 1994]. Comparisons 
between observed and calculated GO are complicated further by 
uncertainties in estimating [Cl y ] and processes that regulate pro- 
duction and loss of CINOj and HG. Stimpfle et al. [1994] 
demonstrate that in situ [CIO] and [N0 2 ] are inversely correlated; 
therefore, errors in the simulation of NO, will be reflected directly 
in the simulation for GO. Significant refinements in atmospheric 


and laboratory observations are required to define the significance 
of the discrepancies for [NO,] and [CIO] illustrated in Fig. 2. 

Concentrations of HG computed using all models are higher 
than measured values by about a factor of 2. Our study demon- 
strates that this discrepancy does not arise from errors in simulat- 
ing HO x . The origin of this discrepancy and its implications for 
understanding the response of 0 3 in the lower stratosphere to 
natural and anthropogenic perturbations remain unclear. 

Effects of NO y Emitted from Aircraft on O 3 

The contributions of HO x , NO x , and halogen catalytic cycles 
to the loss rate of O3 is shown in Fig. 3a-c. Models with hetero- 
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Figure 3. (a-c) The odd oxygen ([0]+[0 3 ]) loss rate, averaged 
over 24 hours, for HO x (dark line; rate limiting step primarily 
H0 2 +0 3 ), halogen (grey line; primarily BrO+CIO, C10+H0 2 , 
B1O+HO2, and CIO+O), NO x (dotted line; primarily N0 2 +0), and 
O+O3 (clashed line) catalytic cycles for the simulations shown in 
Fig. 2. (d) Total loss rate of odd oxygen, averaged over 24 hours, 
for the JPL-Gas (dotted line), JPL-Het (dark line), and Model B 
(grey line) simulations, and the production of odd oxygen from 
photolysis of 0 2 (dashed line), (e) Change in the loss rate of odd 
oxygen (percent) calculated assuming a 20% increase in the con- 
centration of NOj for the JPL-Gas (dotted line), JPL-Het (dark 
line), and Model B (grey line) simuladons, using the observed 
aerosol surface area (Fig. 1). (f) Same as (e), except using esti- 
mates of aerosol surface area for background (non-volcanic) con- 
ditions (-1 x 10" 8 cm" 1 at pressure = 60 mbar). 



2550 


S ALA WITCH ET AL.: DISTRIBUTION OF HYDROGEN, NITROGEN, AND CHLORINE RADICALS 


geneous chemistry on sulfate aerosols indicate larger contributions 
from HO, and halogen cycles and larger total loss rates (Fig. 3d) 
than models based solely on gas phase reactions [Rodriguez et al., 
1991; McElroy et al., 1992]. 

Emission scenarios for the proposed fleet of High Speed Civil 
Transports (HSCT) indicate possible increases in [NO,] of 20 to 
100% at ER-2 flight levels. The gas phase model predicts 
decreases in [Oj] as [NO,] rises due to NO, -related chemistry, in 
contrast to models with heterogeneous chemistry that predict 
increases due to deactivation of the more efficient HO, and halo- 
gen cycles [Weisenstein et al., 1993]. The decrease in the chemi- 
cal removal rate of O3 for an increase in [NO,], for aerosol load- 
ing encountered during SPADE, is -30% larger for Model B than 
for JPL-Het (Fig. 3e), reflecting increased rates for HO, and halo- 
gen cycles in Model B. Loss rates exhibit a more pronounced 
difference to a perturbation of [NO,] for background (non- 
volcanic) aerosol loading due to larger contributions from the NO, 
cycle in each heterogeneous model (Fig. 3f). 

The region of the atmosphere accessible to the ER-2, 16 to 20 
km, is characterized by slow rates for chemical removal of 0 3 . 
Perturbations to [Os] at higher altitudes (20 to 30 km) by HSCT 
aircraft will be sensitive to chemical partitioning of radicals and to 
the transport of NO, emitted from aircraft. Simultaneous meas- 
urements of HO„ NO,, and CIO, radicals and long-lived precur- 
sors at these altitudes are essential for constraining future assess- 
ment models. 

Conclusions 

Models that include heterogeneous hydrolysis of N 2 0 5 on sul- 
fate aerosols provide a good description of observed concentra- 
tions of hydrogen, nitrogen, and chlorine radicals in the lower stra- 
tosphere. Nevertheless, the relative contribution to chemical 
removal of O3 from cycles involving HO, and halogens appears to 
be underestimated by models based on recommended rates and 
cross sections. Consequently, emission of NO, from supersonic 
aircraft may result in larger increases in concentrations of O, in 
the lower stratosphere than predicted by current assessment 
models [e.g., Weisenstein et al., 1993]. 
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